Detected jump-error correcting quantum codes, quantum error designs and quantum 

computation 
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The recently introduced detected-jump correcting quantum codes are capable of stabilizing qubit- 
systems against spontaneous decay processes arising from couplings to statistically independent 
reservoirs. These embedded quantum codes exploit classical information about which qubit has 
emitted spontaneously and correspond to an active error-correcting code embedded in a passive 
error-correcting code. The construction of a family of one detected jump-error correcting quantum 
codes is shown and the optimal redundancy, encoding and recovery as well as general properties of 
detected jump-error correcting quantum codes are discussed. By the use of design theory multiple 
jump-error correcting quantum codes can be constructed. The performance of one jump-error cor- 
recting quantum codes under non-ideal conditions is studied numerically by simulating a quantum 
memory and Grover's algorithm. 



I. INTRODUCTION 

The discovery of powerful, fast quantum algorithms 
||, [| launched new efforts to implement such quan- 
tum algorithms in real physical systems. Quantum al- 
gorithms simultaneously exploit two characteristic fea- 
tures of quantum theory. Namely, the fundamental phe- 
nomenon of quantum interference and the fact that for 
distinguishable quantum systems the dimension of the 
Hilbert space increases exponentially with the number of 
systems. Therefore, to implement a quantum algorithm 
in a real quantum system we must be able to create and 
manipulate arbitrary superpositions of quantum states 
and to preserve quantum coherence during a computa- 
tion. Unfortunately, quantum coherence is very fragile. 
Typically, any coupling to an environment leads to deco- 
herence so that quantum mechanical superpositions are 
rapidly destroyed. 

The urgent need to develop efficient methods to protect 
quantum coherence has led to the study of very general 
classes of quantum error-correcting codes [@, ||, 0, @] • The 
main idea is to restrict the dynamics of a quantum algo- 
rithm to a subspace of the Hilbert space, in which errors 
can be identified uniquely by suitable measurements and 
where the error operators can be inverted by unitary op- 
erations. Typically, this is achieved by an encoding of the 
logical information and by a suitable choice of quantum 
gates. 

For some special cases it is also possible to design a 
passive error-correcting quantum code || ^, [To) . Such a 
passive quantum code relies on a subspace of the Hilbert 
space which is not affected by any errors at all. In this 
situation the unitary recovery operation is the identity 
operation so that an active correction of the errors is not 
necessary. 

In principle, any type of error can be corrected by 



these strategies as long as enough physical qubits are 
available to achieve the required redundancy and one can 
make a large number of control measurements and per- 
form the rapid recovery operations. However, in view of 
current-day experimental possibilities JTT[ it is generally 
difficult to achieve both requirements. Therefore it is de- 
sirable to develop alternative error-correcting strategies 
which possibly correct a restricted class of errors only, but 
which tend to minimize both redundancy and the num- 
ber of recovery operations. Recently, the first steps in 
this direction have been taken by defining a new class of 
one detected jump- error correcting quantum codes which 
are capable of stabilizing distinguishable qubits against 
spontaneous decay processes into statistically indepen- 
dent reservoirs jL2|. These codes are constructed by em- 
bedding an active error-correcting code in a passive code 
space and by exploiting information available on error 
positions. This embedding procedure leads to a signifi- 
cant reduction of redundancy and the number of control 
measurements and recovery operations. 

In this paper the physical principles underlying de- 
tected jump-error correcting quantum codes are explored 
and generalized, motivated by the practical need for 
quantum error-correcting codes which minimize both re- 
dundancy and the number of recovery operations. Based 
on these physical principles an upper bound is established 
on the number of logical states of a general embedded de- 
tected jump-error correcting quantum code. 

From this bound it is apparent that the recently discov- 
ered one detected jump-error correcting quantum codes 
have minimal redundancy. Based on this family of op- 
timal one detected jump- error correcting quantum codes, 
we establish links with the general notions of combina- 
torial design theory fl3| . For this purpose the new con- 
cept of a spontaneous emission error design is introduced. 
This is a powerful tool for constructing multiple detected 
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jump-error correcting quantum codes capable of stabi- 
lizing distinguishable qubits against spontaneous decay 
processes. As an example, we present a new embedded 
three detected jump-error correcting quantum code. 

This paper is organized as follows. In Sec. II basic 
physical aspects concerning the spontaneous emission of 



photons by qubit-systems are summarized. In Sec. Ill 
the physical principles are explored which lead to the 
construction of one detected jump-error correcting quan- 
tum codes. The conditions for general detected jump- 
error correcting quantum codes are given in Sec. IV. The 
links with combinatorial design theory are established in 
Sec. |y|. Finally, in Sec. [vi| numerical examples are pre- 
sented which exhibit basic stability properties of the op- 
timal one detected jump-error correcting quantum codes. 



II. SPONTANEOUS DECAY OF 
DISTINGUISHABLE QUBITS AND PHOTON 
DETECTION 

In this section we summarize basic facts about the dy- 
namical description of a quantum system interacting with 
initially unoccupied modes of the electromagnetic field. 
These considerations are the starting point for the devel- 
opment of optimal strategies of error correction, which 
we pursue in the subsequent sections. 

We consider a model of a quantum computer in which 
N two-level atoms (qubits) interact with external laser 
pulses which synthesize the quantum gates underlying 
a quantum algorithm. These N qubits are assumed to 
be arranged in an array with well defined positions x a 
(a = 1, . . . ,N) (see Fig. [l]). In addition, these N qubits 
are assumed to be distinguishable, which requires that 
their mean nearest neighbor distance is large in com- 
parison with the optical wave lengths involved. Their 
distinguishability guarantees that the dimension of their 
associated Hilbert space is K-n = 2 N and thus scales ex- 
ponentially with the number of qubits. In addition, it is 
assumed that these N qubits couple to the unoccupied 
modes of the electromagnetic field. This coupling causes 
spontaneous decay processes of the qubits from their ex- 
cited states |l) a , to their stable lower lying states |0) a . 
Within the Born, Markov, and the rotating wave approx- 
imations the resulting dynamics of the reduced density 
operator p{t) of this iV-qubit system are described by the 
master equation JTJ, [l^] 



electromagnetic field are described by the Hamiltonian 
H which incorporates the influence of the external laser 
pulses. In addition, we assume that the mean distance 
between the qubits is much larger than the wave lengths 
of the spontaneously emitted radiation. Therefore, to a 
good degree of approximation each qubit a couples to a 
different set of modes of the radiation field so that these 
sets constitute statistically independent reservoirs [p"5| . 
In Eq. ([]) the coupling of qubit a to its reservoir and 
the resulting spontaneous decay process |1) Q — » |0) Q is 
characterized by the Lindblad operator 



L, 



® |0) QQ (1| 



(3) 



where 1/3^ Q denotes the identity on every except the a-th 
qubit, and n a is the associated spontaneous decay rate. 




FIG. 1: Model of a quantum computer consisting of an en- 
semble of distinguishable two level systems (qubits) which are 
well localized at positions x a and which are coupled to the 
vacuum modes of the radiation field. The mean distance d 
between adjacent qubits is assumed to be large in comparison 
with the wave length of the spontaneously emitted radiation, 
i.e. A d. Thus, the spontaneous decay processes do not af- 
fect the distinguishability of the qubits. The qubits are mon- 
itored continuously by photodetectors capable of determining 
the time at which a spontaneous decay process occurred and 
the position of the affected qubit. 



Provided that initially the -/V-qubit system is in a pure 
state, say |~0o), a formal solution of the master Eq. (||) is 
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with the non-hermitian effective Hamiltonian 



H, 



elf 



H 



h N 

a = l 



(2) 



Thereby, the coherent dynamics of the iV-qubit system in 
the absence of the coupling to the vacuum modes of the 



\t;t n a n ,... , txax){t; t n a n , . . . ,hai\ 
with the pure quantum state 
\t; t n a n , . . . , tia\) 

p(t\t n a n , . . . jhaj) 1 / 2 



(4) 
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and with the probabilities 

p(t;t n a n7 ■ ■ ■ ,iiai) (5) 
= (Vole^** 17 ^ . ..Ll n e iH ^-tn)/n . 

It can be shown that each pure state \t; t n a n , . . . , t\ai) 
describes the quantum state of the A-qubit system at 
time t conditioned on the emission of precisely n pho- 
tons at times t\ < ... < t n by qubits ot\, . . . ,a n 
jl4], [H], [l6], [Trl . Thus, each of the pure quantum 
states of Eq. (^|) corresponds to a possible measure- 
ment record . . . ,t n a n ) in an experiment in which 
each qubit is observed continuously by photodetectors. 
In the subsequent discussion it is important to note 
that due to the large separation between the qubits 
ideally this measurement record not only determines 
the spontaneous decay times ti 1 but also the associated 
positions (i = 1, ... ,n) of the qubits which have 
been affected by these decay processes. The measure- 
ment record (t\a.i, . . . , t n a n ) is observed with probability 
p(t;t n a n , . . . , tx<y.\). According to Eq. (Q) the quantum 
state |i; t n a n , . . . , iiai) resulting from a particular mea- 
surement record is determined by two types of effects. 
First, the time evolution between two successive pho- 
ton emission events is characterized by the non-hermitian 
Hamiltonian H e s of Eq. (|2|). Thus, even in the absence 
of any spontaneous photon emission process in a given 
time interval [0, t], in general an arbitrary quantum state 
is modified by the couplings to the vacuum modes of 
the electromagnetic field. Second, immediately after the 
spontaneous emission of a photon by qubit a the quan- 
tum state is modified by the Lindblad operator L a . This 
distinction between the two different types of dynamics 
is crucial for the development of a useful quantum error- 
correcting strategy which minimizes the redundancy as 
well as the number of control measurements and recovery 
operations. 

III. DESIGNING OPTIMAL ONE DETECTED 
JUMP-ERROR CORRECTING QUANTUM 
CODES 

In this section we introduce the main ideas leading 
to the construction of optimal one detected jump-error 
correcting embedded quantum codes which are capable 
of stabilizing systems of distinguishable qubits against 
spontaneous decay processes. These quantum codes ex- 
ploit the physical possibilities which are offered by the 
continuous observation of the distinguishable qubits by 
photodetectors. 

How can we stabilize an A-qubit system, such as the 
one depicted in Fig. ^ against spontaneous decay pro- 
cesses if we can observe the qubits continuously with the 
help of photodetectors? To achieve this two major tasks 
have to be tackled. Namely, we have to correct the mod- 
ifications taking place during successive photon emission 



events. These modifications are described by the effec- 
tive Hamiltonian of Eq. (||). And, we have to invert each 
quantum jump \l) a — > |0) a caused by the spontaneous 
emission of a photon by qubit a. 

In principle, the errors taking place during two suc- 
cessive photon emission events can be corrected by an 
active error-correcting quantum code with a sufficiently 
rapid sequence of control measurements and recovery op- 
erations |l8). However for practical applications such an 
approach is not very attractive. One faces all the prob- 
lems which are related to the implementation of large 
numbers of control measurements and recovery opera- 
tions. It is therefore much more desirable to develop pas- 
sive error-correcting methods for stabilizing the A-qubit 
system against modifications arising from the effective 
Hamiltonian of Eq. (||) . 

The main idea is to determine an appropriate subspace 
D of the Hilbert space of the A distinguishable qubits in 
which the perturbing part of the Hamiltonian acts as a 
multiple of the unit operator, i.e. 

fleffb = H\d - \ct\ D (6) 

where c denotes a c-number. (In our setting this number 
is real- valued and positive.) If the (unperturbed) dy- 
namics characterized by the Hamiltonian H do not take 
an initial pure quantum state outside this decoherence 
free subspace D, this condition implies that, provided no 
photon is emitted in the time interval [0, t], the quantum 
state at time t is identical with the unperturbed state 
(compare with Eq. (|4|)). Thus, if one can find such a suf- 
ficiently high dimensional decoherence free subspace, the 
dynamics taking place between successive spontaneous 
photon emission events are stabilized perfectly without 
the need for control measurements and recovery opera- 
tions. In practice it is desirable to choose the dimension 
of the decoherence free subspace to be as large as possi- 
ble. 

An important special case occurs when all the A qubits 
have identical spontaneous decay rates, 1.6. — t\j . in 
this situation it follows that 

N N 

£4L a = K£|l) Qa <l| (7) 

OL — l OL — 1 

and any subspace formed by basis states involving an 
equal number, say w < A, of excited qubits is a deco- 
herence free subspace. For a given number of A qubits 
the dimension Kd of such a decoherence free subspace 
is given by Kjj = ( N ) which is maximal if w = \_N/2\. 
(\_x\ denotes the largest integer smaller or equal to x.) 

In general, the first spontaneous emission of a pho- 
ton will affect the quantum state of the A-qubit system 
in an irreversible way. According to Eq. (^) the spon- 
taneous emission of a photon by qubit a, for example, 
is described by the application of the Lindblad operator 
L a which induces a quantum jump |1) Q — > |0) Q . This 
Lindblad operator is not invertible over the decoherence 
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free subspace D so that this quantum jump cannot be 
corrected. In order to correct for this quantum jump ac- 
tively we have to restrict the dynamics to a still smaller 
subspace C C D in which a unitary operator, say U a , 
can be found having the property 

U a L a \ c = l| c . (8) 

Therefore, if we still want to take advantage of passive er- 
ror correction between successive photon emission events 
we have to construct an active error-correcting quantum 
code within the relevant decoherence free subspace. 

We now construct a one detected jump-error correcting 
embedded quantum code in the special case of identical 
spontaneous decay rates considered above. According 
to the criterion given in the orthogonal basis states 
{|cj) : i = 1, . . . , K} of a subspace C constitute an active 
error-correcting quantum code with respect to the set of 
error operators {L a : a = 1, . . . , A} if and only if 

{a\LlLp\ Cj ) = %A Q/3 (9) 

for all possible values of i,j and a, f3. Eq. (||) states the 
necessary and sufficient conditions for the existence of 
unitary recovery operations which fulfill Eq. (^]) for the 
error operators {L a : a — 1,... , A}. In the physical 
setting this criterion states that: (i) different orthogo- 
nal quantum states remain orthogonal under the action 
of error operators; and (ii) all basis states are affected 
by a given pair of errors L a and Lp in a similar way. 
The latter condition necessarily implies that the scalar 
products between the states L a \ci) and Lp\d) are state- 
independent. It is plausible that a larger set of error 
operators leads to a more restrictive set of conditions of 
the type of Eq. (^|). Furthermore, we also expect that 
more restrictive conditions lead to a higher redundancy 
of an active quantum code. 

As an example, consider the situation where contin- 
uous observation of the A-qubit system by photodetec- 
tors does not reveal which qubit a has emitted the regis- 
tered photon. This implies that the error operators which 
could induce a spontaneous decay process are in the set 
{L a : a = 1, . . . , A}. It has been shown by Plenio et 
al. [|l9| that when the error positions are unknown, eight 
physical qubits are needed to encode two orthogonal log- 
ical states by an embedded quantum code. This should 
be compared with the optimal active one-error correcting 
code using five qubits |2(J . 

Thus, the advantage offered by using an embedded 
quantum code, capable of passively correcting errors be- 
tween successive photon emission events leads to a sig- 
nificant increase of redundancy in comparison to purely 
active methods. However, this disadvantage can be over- 
come if besides knowing the error time, information 
about the error position a is also available. In principle, 
this information can be obtained from continuous obser- 
vation of the A-qubit system by photodetectors as long 
as the mean distance between adjacent qubits is large in 
comparison with the wave length of the spontaneously 



emitted radiation. For this purpose it is important that 
each photon which is emitted by one of the qubits can be 
detected. 

How can we construct a one detected jump-error cor- 
recting embedded quantum code which exploits infor- 
mation about the error position in an optimal way so 
that its redundancy is minimized? Let us concentrate 
again on our previously introduced example of identi- 
cal spontaneous decay rates. In this setting we have 
a decoherence free subspace which involves w excited 
qubits. This stabilizes the dynamics between successive 
photon emission events passively. For example, in the 
simple case of A = 4 and w — 2, the orthogonal ba- 
sis states |cj) of the decoherence free subspace are given 
by { 1 1100), |0011), |1010>, |0101), |1001), |0110)}. Within 
this six-dimensional decoherence free subspace a possible 
active quantum code, capable of correcting one quantum 
jump at a time at a known position a, is determined by 
the following three orthonormal code words: 

|c ) = i[|1100) + e^|0011)], 

|d) = i[|01010) + e^|1010)], 

|c 2 ) = i=[|1001) + e^|0110)]. (10) 

These code words are formed by all possible different 
pairs of basis states of the decoherence free subspace 
which can be constructed by interchanging states |0) and 
|1). These complementary pairings are also character- 
ized by equally weighted probability amplitudes which 
involve an arbitrary phase ip. It can be easily checked 
that this embedded quantum code fulfills the criterion of 
Eq. (||) for any of the Lindblad operators L a of Eq. (||) 
separately, but not for two different Lindblad operators 
with a ^ (3. Thus, provided that the error time and the 
error position are known this embedded quantum code 
can correct one quantum jump at a time. In addition, 
errors between successive quantum jumps are corrected 
passively. In this way the number of control measure- 
ments and recovery operations is reduced and so is the 
redundancy. This is apparent by comparing our quan- 
tum code with the code proposed by Plenio et al. [ p^[ , 
which requires eight physical qubits for the encoding of 
two orthogonal quantum states. 

The above construction based on complementary pair- 
ings can be generalized to an even number A of physical 
qubits; giving an infinite family of one detected jump- 
error correcting embedded quantum codes. The con- 
struction is the following. We start from the highest di- 
mensional decoherence free subspace which involves all 
quantum states in which half of the qubits are excited to 
states |1) Q . Subsequently, we construct an active quan- 
tum code within this decoherence free subspace by com- 
plementary pairings. The dimension Kj of the resulting 



5 



embedded quantum code is given by 



Kj = (1/2) 



N 
N/2 



N - 1 
N/2 - 1 



(11) 



It is a remarkable fact that, provided one wants to correct 
errors passively between successive quantum jumps, it is 
not possible for an even number of physical qubits N to 
reduce the redundancy of such an embedded jump code 
(see Sec. and Ref. ||). Therefore, our family of 
embedded quantum codes has the desirable feature that 
it minimizes redundancy and that it requires only a small 
number of recovery operations. For instance, in a given 
time interval, say [0, i], the mean number n of required 
recovery operations is of the order ofnw 2t/(Nn). 

To stabilize any quantum algorithm against sponta- 
neous decay processes using an embedded one detected 
jump-error correcting quantum code three requirements 
have to be met. First, one has to be able to register the 
time t and position a of each spontaneous decay event 
which takes place during the performance of the quan- 
tum algorithm. As indicated schematically in Fig. |l] this 
can be achieved by continuous observation of the iV-qubit 
system with photodetectors. In principle, an identifica- 
tion of the perturbed qubit a is possible provided the 
mean nearest neighbor spacing of the qubits is large in 
comparison with the wave lengths of the radiation emit- 
ted spontaneously. However, in practice the error posi- 
tion a might not be determined so easily due to imperfect 
detection efficiencies of the photodetectors. Therefore, in 
actual applications shelving techniques [^2], might be 
useful which amplify each spontaneously emitted photon 
to such an extent that it can be detected with an effi- 
ciency arbitrarily close to unity. 

Second, we have to ensure that each spontaneous de- 
cay event is corrected immediately by application of an 
appropriated unitary transformation which inverts the 
effect of the Lindblad operator L a . In practice, this in- 
version has to be performed on a time scale which is 
small in comparison with the natural time scale of the 
quantum algorithm and with the mean spontaneous de- 
cay time. Third, one has to ensure that the sequence of 
quantum gates which constitute the quantum algorithm 
does not leave the code space at any time. This can be 
done by encoding the logical information within the code 
space, and to develop a universal set of quantum gates 
which leaves this code space invariant. 

Ideally these quantum gates are implemented by suit- 
able Hamiltonians. This ensures that the code space 
is left invariant even during the application of one of 
these universal quantum gates. Such universal sets of 
Hamiltonian-induced quantum gates have already been 
developed for dccohercncc free subspaces D of the kind 
discussed above || [|, [To). But in general, unitary gates 
based on swapping Hamiltonians need not be universal on 
the embedded quantum code, or the swapping Hamilto- 
nians do not leave the embedded quantum code invariant. 
The solution of this intricate and yet unsolved problem 
is beyond the scope of the present work. However, some 



pre liminary results have already been obtained recently 

So far we have shown that any Lindblad operator L a 
of the form of Eq. (||) can be inverted by our one de- 
tected jump-error correcting quantum codes. We provide 
an example of a unitary transformation which achieves 
this inversion in the case of the one detected jump-error 
correcting quantum code involving four physical qubits. 
A possible sequence of quantum gates capable of invert- 
ing a spontaneous decay process affecting qubit a = 2, 
for example, is depicted in Fig. |2|. This example demon- 
strates the basic fact that it is indeed possible to perform 
a unitary inversion of the Lindblad operator Li provided 
Eq. (p|) is fulfilled for a = (3 = 2. 
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Classical control of recovery 



FIG. 2: Gates to invert a quantum jump (indicated by a bolt) 
of a physical qubit. The four physical qubits involved form a 
(4, 3, 1) 2 code. 



IV. GENERAL DETECTED JUMP-ERROR 
CORRECTING QUANTUM CODES 

To define general detected jump-error correcting quan- 
tum codes, we introduce some notation. For a set of 
positions E = {a\,a2, ■ ■ ■ ,cv„}, we denote by Je the 
operator 



Je '■= J{ 



The associated error times t\, . . . ,t n are no longer men- 
tioned explicitly, but it is understood that they are 
known. Note that the operators L ai , . . . , L an commute, 
because the cti are pairwise different. Since by Eq. (|^), 
the errors which involve two equal indices, say ot\ =012, 
cannot occur. 

As discussed in Sec. Ill, for H = all the states which 



are superpositions of states with a constant number of ex- 
cited qubits are common eigenstates of the non-unitary 
effective time evolution (^) between quantum jumps. A 
subspace C of such a decoherence free subspace with or- 
thonormal basis {|ci),... , \ck)} is called a d-detected 
jump-error correcting quantum code, and is denoted by 
C = (N, K, d) w if the following condition holds for sets 
E of jump positions with at most d elements and for all 
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basis states \q) and \a 



(c i \J E J E \c j )=6 ij \{E), 



(12) 



The notation C — (N,K,d) w is motivated by classical 
coding theory. Similarly, the notation d-JC(N, w, K) is 
motivated by notations from design theory (see |fL2[ ), 

The validity of this statement follows from the general 
conditions on quantum error-correcting codes (cf. ||). 
Since we know on which positions the jump operator J E 
acts, only products of the form 3 e Je have to be consid- 
ered. 

There is a natural connection with combinatorics. For 
a basis state \x\ . . . xn) of N qubits, the positions which 
are in state |1) define a subset of {1,... ,N}. Fur- 
thermore, a collection of such subsets corresponds to an 
equally weighted superposition of basis states. 

Let SW, . . . ,B^ be K disjoint sets of subsets X C 
{1, . . . ,N}, where each subset X contains w elements. 
Identifying the set X and the binary word x where Xi = 1 
if i £ X and X{ = otherwise, we define the states 



1 



Ci) = 



XGBW 



(13) 



Where \B"'\ denotes the number of elements (cardinal- 
ity) of the set $W. These orthonormal states span 
a d-detected jump-error correcting quantum code C = 
(N, K, d) w provided that for all sets of jump positions 
E with no more than d elements and all sets B^ the 
following condition holds: 



\X £ gM : E C X}\ 



= KE)- 



(14) 



We note that the disjointness of the sets B"' implies con- 
dition Jl2] ) for i ^ j. Rewriting the operator J e Je as 
J e Je = J2xde \ x )( x \ shows that for i = j and the 
states ( |T3| ) the expectation value (ci\J E JE\<k) equals the 
expression in (14). 



V. BLOCK DESIGNS AND JUMP CODES 

In this section we will show that ci-detected jump-error 
correcting quantum codes are naturally connected with 
(i-designs. These are combinatorial structures which have 
been extensively studied for many decades, cf. |jl3). To 
denote this class of combinatorial structures we intro- 
duce some notations using the language of finite inci- 
dence structures. Let V be a set of N elements, called 
points, say V — {1, ... ,N} where N > is an integer 
and |V| = N. For < w < N the class of w-subsets of 
V containing w elements will be denoted by ( J . In a 
suggestive way its cardinality, i.e. 



V 



w 



N 



is just the binomial coefficient ( ) . An incidence struc- 
ture in V is specified by a distinguished class B of subsets 
of V. The elements B of B are called blocks (or some- 
times lines) of the incidence structure. If B C ( J , we 
say that B has constant block size w. As an example, 
any undirected graph is an incidence structure of block 
size two, if we choose V as the set of points of the graph 
and the points which are directly connected by an edge 
as blocks. 

For any point a £ V the class B a := {B G B : a £ B] 
denotes the class of blocks containing the point a (or: 
"the lines through a"). If \B a \ = r is constant for all 
a £ V the incidence structure B is called regular. For 
a graph, \B a \ is the degree d(a) of the vertex a, i.e., 
the number of edges on which a lies. The incidence 
structure B — as well as the graph itself — is called reg- 
ular, if d(a) is constant. If there exists a constant 
such that for all E = {ai,... ,ctd} £ (v) the class 
Be ■= (((Bai)a 2 ) ■ ■ - )a d of all blocks containing the set 
d of points has size \Be\ — ^d, the incidence structure 
is called d-regular. A d- design is a d-regular incidence 
structure with constant block size w. It is denoted by d- 
(N, w, Ad)-design or as S\ d (d,w,N) p3| . Graphs which 
correspond to 2-designs are depicted in Fig. ||. 





Si (2, 2, 4) 



Si (2, 3, 7) 



FIG. 3: Examples of 2-designs. 

The preceding discussion leads to the notion of d- 
spontaneous emission error designs which we denote by 
d-SEED(N,w,K) The essential property of these 

combinatorial objects is the local multiplicity X(E) of a 
subset E of {1, ... ,N} containing at most d elements 
which is defined by 



JBf| 
\BW\ 



(15) 



where B^\... are disjoint subsets of (^). Any 

d-SEED(7V, w, K) produces a d-detected jump-error cor- 
recting quantum code C = (N, K, d) w using the encoding 
defined in Eq. @. 

We conclude this section by constructing of a three- 
jump correcting code. The permutation group G — 
((1 2) (3 4), (1 4)(2 3), (5 6)(7 8), (5 8)(6 7), (1 2 3)(5 6 7)) of 
order 48 acts on the 4-element subsets of {1,2,... ,8}. 
The orbits under G of the sets {1,2,5,6}, {1,3,5,6}, 
and {1,4,5,6} are mutually disjoint. Direct calcula- 
tion shows that they fulfill the local multiplicity condi- 
tion @. Hence the sets B^ := {1,2,5,6} G , B^ := 
{1,3,5,6} G , and B^ := {1,4,5,6} G define an (8,3,3) 4 
jump code. The corresponding (not normalized) basis 
states are given by: 
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|ci) = 100110011} + 100111100} + |01010101} + 101011010} 

+ 101100110} + 101101001} + |10010110} + 110011001} 
+ |10100101} + |10101010} + 111000011} + 111001100} 

|c 2 } = 100110110} + 100111001} + 101010011} + 101011100} 

+ 101100101} + 101101010} + |10010101} + 110011010} 
+ |10100011} + 110101100} + 111000110} + 111001001} 

|c 3 } = 100110101} + 100111010} + 101010110} + 101011001} 

+ 101100011} + 101101100} + |10010011} + 110011100} 
+ |10100110} + |10101001} + |11000101} + 111001010} 

Further examples of d-SEEDs are discussed in EJ . In 
that paper, there are also general bounds on the param- 
eters of jump codes derived. In particular, the dimension 
if of a d-detected jump-error correcting (N, K, d) w code 
is bounded above by 



K < min{ 



N ■ 



w 



N-d 



w 



}■ 



(16) 



For completeness, we repeat the main ideas of the proof. 

The dimension of the space spanned by the basis states 
of N qubits where w of them are in the excited state is 
( U) ). This implies the bound K < A jump on 

j positions reduces the number of excitations to w — j. 
After the jump, the j positions where the jump occurred 
are zero. There are (^Zj) such basis states. A jump must 

not reduce the dimension of the code, hence K < ( N ~ : ') . 

5 — \w—j/ 

For < j < d possible quantum jumps the lowest upper 
bound is achieved for j = d, as 



N 



w 



N(N - 1) . . . (N - d + 1) (N - d 



w(w — 1) . . . (w — d + 1) 



-d 



To obtain the second upper bound in Eq. (|l6|), we note 
that starting with a jump code C — (N, K,d) w , applying 
u x to all qubits yields a jump code C — (N,K,d)jy^ w . 
Note that, a x interchanges ground and excited state, 
hence the code C lies in the decoherence free subspace 
with N — w excitations. Also, the linear space spanned 
by the operators J\Je for all subsets E with no more 
than d elements is invariant under conjugation by a x 
on all qubits. This holds as a x L^< a L a <7 x = |0) aQ (0| = 
1 — IDn-rvdl- Thus, for the code C we obtain the bound 



K < 



' N-d 
.N—w- 



,) = ( ) . If there is no restriction on the 
number w of excited states, choosing w — \_N/2\ maxi- 
mizes the upper bound of Eq. (jig), i.e. 



K<( N ~ d 
~ \\N 2\-d 



(17) 



As mentioned above, the upper bound of Eq. (^) is 
achieved for d — 1 and for an even number N of qubits. 
A table of lower bounds (obtained by constructions from 
rf-SEEDs) and upper bounds for small values of JV and d 
is provided in plL 



VI. STABILITY PROPERTIES OF THE ONE 
DETECTED JUMP-ERROR CORRECTING 
QUANTUM CODES 

The detected jump-error correcting quantum codes 
constructed in the previous sections can stabilize quan- 
tum algorithms provided three conditions are fulfilled. 
First, the decay rates of the qubits are equal. Second, 
the time and position of each quantum jump is detected 
with hundred percent efficiency. Third, the appropriate 
unitary recovery operations are applied perfectly and in- 
stantaneously immediately following the detection of an 
error. Experimentally these requirements can only be ap- 
proximated. Therefore, the natural question arises how 
non-ideal conditions influence the robustness of the em- 
bedded error-correcting quantum codes. 

In this section several types of imperfections affecting 
the ideal performance of the one detected jump-error cor- 
recting embedded quantum codes of Sec. Ill are studied 
numerically. For this purpose we investigate the stabi- 
lization of a quantum memory and of a simple Hamilto- 
nian dynamics against spontaneous decay processes. The 
effective two-level Rabi Hamiltonian considered, i.e. 



H = iMl (\x )(v\-\v)(x \), 



(18) 



can be viewed as modeling the quantum dynamics of 
the ideal Grover search algorithm in the limit of a 
large number of qubits. Thereby, the Rabi frequency 
Q can be related to the characteristic time r op required 
for performing an oracle-operation and to the number 
of qubits N according to SI = 2/(t op VN) Here, 
\v) denotes an equally weighted superposition of all (or- 
thonormal) code words which may be used as an initial 
state in Grover's quantum search algorithm. The final 
state we are searching for is denoted \xo). For this choice 
of \v) the states \v) and \xq) are not orthogonal because 
(v\xq) — 1/yN. However, if the number of qubits N be- 
comes large their overlap tends to zero. According to the 
Hamiltonian Eq. (|l^) and consistent with Grover's quan- 
tum algorithm [2j, after the interaction time t — n/ (20), 
the initial state |i>) is transformed to the final state \xq). 



A. Imperfect detection of error positions 

Let us first of all consider situation where the jump 
position can be detected with a given non-zero error rate 
only. Such an imperfection might occur if, for example, 
a photon emitted by a particular trapped ion is detected 
by the photodetector associated with a different ion (see 
with Fig. |). 

The probability to detect the emitted photon at the 
correct position is denoted by p. The probability that 
an emitted photon is detected falsely by the next nearest 
neighbor is given by q x p. Analogously, the probability 
of detecting the photon by the n-th nearest neighbor is 
q" x p with the normalization condition px (1 -\- q + q 2 + 
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FIG. 4: Model for a 4-qubit register where jumps are de- 
tected with a given error parameter q. 



The influence of this type of imperfection on a quan- 
tum memory, i.e. H = in Eq.(jl]), is depicted in Fig. |^. 
A state IV'(O)) of the jump code C — (4, 3, 1)2 is propa- 
gated according to Eq. (Q) with a Monte-Carlo simulation 
of the quantum-trajectories [p"7| . If a jump is detected, 
the appropriate recovery operation is applied. In the case 
of a correct detection, the quantum state of the memory 
is recovered perfectly. In the case of a false detection, 
the quantum state of the memory leaves the code space. 
Therefore, in this simulation the full Hilbert space of four 
physical qubits has to be taken into account. As a mea- 
sure of fidelity the squared absolute value of the overlap 
between the state IVK 7 ")) after a time r = 7t/(2k) and the 
initial state \ip(0)} of the memory is plotted as a function 
of the error parameter q. 
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ratory, it is of interest to investigate what happens if this 
condition of equal decay rates is violated. In this latter 
case our code does not correct errors between successive 
quantum jumps passively. 

For this purpose let us consider the Rabi-hamiltonian 
of Eq. ( fl8| ) which describes the ideal quantum dynamics. 
In addition, we assume that the decay rates of the physi- 
cal qubits are selected randomly according to a Gaussian 
distribution whose mean value k is equal to the charac- 
teristic Rabi frequency Q. 

To study the resulting time evolution we choose the 
(6, 10, 1)3 code based on complementary pairings (see 
Sec. In Fig. || the fidelity of the quantum state 

p(r) is depicted as a function of the variance An of the 
Gaussian distribution. The fidelity is defined as the over- 
lap (xq\p(t)\xq) between the (mixed) system state p(r) 
at time r = 7r/(2f2) and the desired state |xo)(a;o| which 
would result from the ideal dynamics at this particular 
interaction time. In this numerical simulation the mas- 
ter equation ([!]) was integrated up to time r = 7r/(2f2), 
whereas each jump operator L a was replaced by a se- 
quence consisting of L a and an immediately applied uni- 
tary recovery operation U a (see Eq.(|])). In this simu- 
lation it was assumed that the recovery operations are 
performed perfectly. It is apparent that the code sta- 
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FIG. 6: Dependence of the fidelity of an effective Rabi oscilla- 
tion at time r = 7r/(2fi) on the variance Ak of the decay rates: 
The mean value of the decay rates is k = £1, the variances are 
chosen from Gaussian distributions with Ak = 0.1 ... 1.0 X O. 
The diamonds show the simulation using the one detected 
jump-error correcting quantum code (6, 10, 1)3, the stars show 
the corresponding results for six qubits without encoding. 



B. Unequal spontaneous decay rates 

For n a = k the orthonormal basis states |c,) of the 
code C C D are degenerate eigenvectors of the opera- 
tor Y] L^ a L a appearing in the effective Hamiltonian of 
Eq. (|2j). This property ensures that these states form 
a passive code for the effective time evolution between 
successive quantum jumps. The existence of such degen- 
erate eigenstates of the operator ^ Q L^La relies on the 
assumption that the decay rates of all qubits are equal. 
Although this physical situation can be realized in a labo- 



bilizes the quantum dynamics successfully, despite the 
fact that the (6, 10, 1)3 code is not a perfect one detected 
jump-error correcting quantum code for this situation. 

C. Time delay between error detection and 
recovery operation 

Immediately after the detection of a spontaneous emis- 
sion event the qubits are described by a quantum state 
belonging to a subspace involving one excited qubit less 
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than the original code space. This subspace also consti- 
tutes a passive error-correcting code. Therefore, a time 
delay between the detection and the application of a re- 
covery operation does not lead to an additional error 
caused by the effective time evolution H c g provided the 
ideal quantum dynamics characterized by H is not af- 
fected. Nevertheless, this time delay must be short in 
comparison with the mean time between two successive 
spontaneous emission events. Otherwise, a second spon- 
taneous emission may map the state of the system onto 
another subspace, from which a recovery is no longer pos- 
sible. Fig. [7] demonstrates that, as long as the delay be- 
tween detection and correction is not too large compared 
with the mean decay time, error correction is still possi- 
ble. 
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FIG. 7: Fidelity |(^)(r)|a;o}| 2 for an effective Rabi oscillation 
at time r = 7r/(2f2) with a delay time t between detection 
and recovery: The decay rate is k = 0.5 Q. 



D. Dead times of the photodetectors 

Another important condition for correct implementa- 
tion of a detected jump-error correcting quantum code 
is the ability to observe the environment of each qubit 
continuously. However, immediately after the detection 
of a spontaneous emission event, typically the detector is 
not able to respond to another photon. During the latent 
response time of a photodetector a second spontaneous 
emission event can take place which may destroy quan- 
tum coherence. In Fig. the dependence of the fidelity 



T )\ x o)\ on the response time of a photodetector is 



depicted for various decay rates n. 

It is apparent from Fig. ||that the detected jump-error 



correcting quantum code can stabilize an algorithm as 
long as the response time of the photodetectors is small 
in comparison with the average time between successive 
spontaneous emission events. 
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FIG. 8: Dependence of the fidelity | (^(t) |a;o) | 2 at time 
t — 7r/(2fi) on the dead time t of a photodetector. The 
ideal dynamics are governed by an effective Rabi oscillation 
with Rabi frequency f2. For the encoding the one detected 
jump-error correcting quantum code (6, 10, 1)3 is used. The 
fidelity is depicted for four different values of H: k — 0.25 
(black diamonds), k = 0.5 SI (triangles), 7c = 0.75 (white 
diamonds), k = 1.0 Q (crosses). 



VII. CONCLUSIONS 

We have studied quantum error-correcting codes that 
exploit additional information about the locations of the 
errors. This information is obtained by continuously 
monitoring the system. Errors caused by the resulting 
non-unitary dynamics are corrected passively by embed- 
ding the error-correcting code in a decoherence free sub- 
space. To construct such codes we have established con- 
nections to design theory. The numerical simulations pre- 
sented demonstrate that the jump codes discussed can 
stabilize quantum systems even in cases of imperfect de- 
tections and recovery operation. 
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